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Abstract 

We develop a Lagrangian approach for constructing a symplectic structure for 
singular systems. It gives a simple and unified framework for understanding 
the origin of the pathologies that appear in the Dirac-Bergmann formalism, 
and offers a more general approach for a symplectic formalism, even when 
there is no Hamiltonian in a canonical sense. We can thus overcome the 
usual limitations of the canonical quantization, and perform an algebraically 
consistent quantization for a more general set of Lagrangian systems. 
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I. INTRODUCTION 



There is a canonical recipe for accomplishing the conventional quantization of a dynam- 
ical system. From a given Lagrangian we use a Legendre transformation to rewrite the 
velocities in terms of the momenta and to determine the Hamiltonian, and next from the 
Poisson brackets between coordinates and momenta we define the quantum commutation 
relations. This procedure works reasonably well for regular Lagrangians, but it becomes 
troublesome for the singular ones, which are characterized by the non-invertibility of their 
Hessian matrix. This not only prevents us from determining all the accelerations in terms 
of the coordinates and their velocities, but it also confronts us with a singular Legendre 
transformation when we try to construct the canonical theory. These difficulties reflect the 
use of an excess of coordinates in the original description, because the true physical degrees 
of freedom are contained in a subspace of the configuration space. This subspace is defined 
by a set of Lagrangian constraints, and its complement is a set of physically equivalent 
configurations or a physically inaccessible region. 

Some decades ago, Dirac and Bergmann and colaborators made a very relevant contribu- 
tion to the subject of singular Lagrangians, and developed a systematic way of constructing 
a canonical formalism for these systems The Dirac method has played a crucial role 

in the quantization of many theories of significant physical relevance, as is the case of gauge 
theories and their gauge fixing procedure. However, there exists some pathological exam- 
ples where the Dirac method seems to fail apparently because the constraint analysis 
on which it is based is not exhaustive enough In these cases the Dirac construction 

of the phase space leads to some inconsistencies in relation to the original Euler-Lagrange 
equations. The meaning of the primary/ secondary classification of the constraints is not 
clear 0,^, and the connection between the constraints and the symmetries of the system is 
rather vague. This has been manifested by a large body of work speculating on the use of 
full or partially extended Dirac Hamiltonians, following more or less ad hoc criteria ||5|-0. 
Besides this, for Lagrangians which are not bilinear functions of the velocities, a one to 
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one relation between both velocities and momenta might not exist. This introduces some 
additional problems with the use of the Legendre transformation and spoils the construction 
of a canonical classical formalism and of a consistent canonical quantization PJTO|]. If we 
consider the Lagrangian description for a dynamical system as the basic one, these prob- 
lems are related to the fact that the Dirac approach gives a set of constraints which are not 
unambiguously defined and with no direct relation with the Lagrangian constraints. 

The Dirac approach is an extension of the most usual method of constructing a canonical 
formalism for regular systems, based on a Legendre transformation. But this is not the only 
approach to achieve this aim. One possiblity, based on the theory of transformations for 
the Lagrangian formalism, has been explored in Ref. |Tl|]. Another possibility is to visualize 



the canonical formalism for regular systems as a very special first order Lagrangian theory. 
This first order theory is constructed from the original one by redefining the velocities as 
new auxiliary variables, which leads to a new configuration space of double dimension. This 
is the perspective that we exploit in the present work, which allows us to obtain a canonical 
formalism maintaining a clear relationship with the original Lagrangian, and avoiding the 
ambiguities involved in the usual construction of the phase space. 

Here we develop, in a heurisitic language and using local coordinates, such an alternative 
procedure. It is based on a redefinition of the original velocities as new variables, leading 
to a Lagrangian with a linear dependence on the new velocities in a configuration space 
of double dimension. Essentially, it is an extended version for singular systems of the 
approach discussed by Lanczos [|1^] to accomplish the transition from the Lagrangian to the 



Hamiltonian form, which is based solely on the method of Lagrangian multipliers. In this 
context we construct a symplectic structure for singular Lagrangians. The use of Lagrange 
multipliers allows us to overcome the sometimes not well defined Legendre transformation. 

A few years ago, Faddeev and Jackiw ||T^ proposed a method, based on the classical 
geometrical approach to dynamical systems, in order to identify cyclic variables and the 
associated constraint functions. They then proceeded to use them to eliminate superfluous 
degrees of freedom. However, it sometimes happens that such an elimination is not conve- 



nient or even possible. To handle these situations, these constraints have been incorporated 
into the kinetic terms by using velocities as Lagrange multipliers [0]. This procedure has 
been shown to be equivalent to the Dirac method when the latter is applicable [jT5[, and 
offers a clear description of the gauge symmetries |]T6| . Here we apply the schema presented 
in Ref. |T^-|IB| to develop a transparent analysis of the structure of singular Lagrangians 
and their constraints in a generalized symplectic context. In particular, this analysis sheds 
light on the Dirac conjecture, i.e., on the relationship of first class secondary constraints 
with gauge symmetries [|l^ . 

A very important feature of our approach is that the resulting symplectic formalism 
guarantees the soundness of the original Lagrangian equations of motion. We will deal here 
with finite-dimensional systems, but the main results can be extended without difficulty to 
field theory. 

The development of the article is as follows. In Section II we discuss how to construct 
from the original Lagrangian a first order one that leads to a presymplectic structure. If 
the original Lagrangian is regular, at this stage we already have the Poisson brackets and 
the symplectic formalism. In the three following sections we follow the analysis for singular 
Lagrangians and discuss the appearance of constraints. Here we use the derivative Lagrange 
multipliers to incorporate the constraints to the first order Lagrangian and study the different 
cases. The resulting schema allows us to give a complete characterization of the symmetries 
in relation to the constraints, and to construct a well defined symplectic formalism. To 
illustrate our approach. Section VI contains two examples, already discussed in the literature 
as presenting a pathological character from the Dirac point of view. Here we show how our 
approach allows us to treat them. Finally, the last section is devoted to some comments and 
conclusions. 
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II. THE PRIMARY FIRST ORDER LAGRANGIAN AND ITS CONSTRAINTS 



Let us consider a dynamical system described by a Lagrangian L{qi,qi), of arbitrary 
order in the time derivatives (ji, i = 1, ...,n. To develop a symplectic form for this system, 
we start by transforming this Lagrangian to an adequate first order one, by extending the 
schema of Lanczos |]T2| to singular Lagrangians. The procedure consists in substituting in 
the Lagrangian the velocities Qi by n auxiliary variables ctj to define 

L{qi,ai) = L{qi,q,)\q^=a^ , (1) 

and in constructing a new Lagrangian L^^^{q,q,iT,a), by adding to L{q,a) the constraints 
ai = (ji through a set of Lagrange multipliers TTf. 

^(0) TT, a) = n,{qi - ai) + L{q, a) . (2) 

This new Lagrangian leads to exactly the same equations of motion for the variables 
than the original L{qi, q^), as can be easily verified by eliminating the cyclic variables (vTj, a,) 



T7| . The structure of L^^^ becomes more meaningful when isolating the velocity- independent 



terms: 

-^(0)(9>9>7r,a) = TTiQi - i7(g,7r,Q;) . (3) 

The function of the extended configuration space 

H{q,n,a) = aiHi - L{q,a) (4) 

is a constant of motion for the Lagrangian L^^^{q,q,7T,a). The are cyclic variables, and 
their equations of motion give rise to a set of primary constraints: 

These constraints resemble the definition of the canonical conjugate variables in the usual 
Legendre transformation, but here they are defined in an extended configuration space, and 
simply characterize its enlargement. Expression @ also strongly resembles the usual relation 



between the Lagrangian and the Hamiltonian function but here, due to the constraints, the 
equations of motion cannot be cast into a canonical form and it is not possible to define a 
kind of generalized Poisson brackets. 

On the subspace defined by (|^) the function H is a -independent. If the constraint 
= allows us to algebraically explicitate a^, it can be used to eliminate this variable 
from the Lagrangian L^^^ [Q, and thus to reduce the dimension of the extended configuration 
space. By doing this we show in Appendices A and B how the canonical and the Dirac 
Hamiltonian appear. 

In the following, to make clear the role of the Lagrangian constraints, we will maintain 
the whole set of auxiliary variables. To promote the constraints to Euler-Lagrange equations 
they must be incorporated to the Lagrangian L^^^ {q, q, ir, a) by means of Lagrange multipli- 
ers. We characterize the constraints by their time derivatives, instead of the constraints 

' (1) 

themselves |]l4| , p!5| . This is easily implemented by introducing velocities as Lagrange multi- 
phers, A*^^ , in place of the standard ones. This procedure has the advantage that only the 
velocity-dependent term is modified, while the velocity-independent part remains the same. 
Thus, from here on the resulting function H{q, tt, a) is completely defined, and the insertion 
of the new constraints will only modify the coefficients of the velocities in the equations of 
motion. In such a way we obtain a new Lagrangian L^^j (g, q, ir, a, A^^j): 

L^,^ {q, q,TT,a, A^^ ) = yr^^i + fi*^^ AJ^^ - H{q, it, a) . (6) 

For the sake of brevity, let us introduce the compact notation Q^^^ = (g; vr; a; A^^^), 
A'^^ = .^^ , where the range of the A index is {l...n; such that L^^^ may 

be written as 

L^,,{QJ = AlX^-H{QJ, (7) 

with H{Q^^^) = H{q,iT,a) given by Eq. (H). This Lagrangian leads to the first order 
equations of motion: 

^(1) 



where F, , = 



(1) 



(1) dQ^^^ dQ^^ 



y- has the exphcit form 



(1) 



/ 













R, 


-K 








(1) 
(1) 



(9) 



The matrix / is the symplectic matrix in the {q, tt) space 

I 



where / is the unit matrix, and the remaining components are 

a 
an 



M, 



(1) 



(1) 



dq 



(1) 



da 



(10) 



9^L d'^L(q,a) 

We can see that the matrix R,^,, whose elements are Rl\ = ^ = — — — , is the Hessian 

^ ' oai oaiOaj 

of L{q, q) evaluated in q = a. 

To make apparent when F^^^ is regular or singular we can perform the following manip- 
ulation. Introducing the idempotent matrix 



^ I f-'M^^^ ^ 



U, 



(1) 



/ 
-/ 



(11) 



we construct 



W F U 



/ 
-R 



(1) 



(12) 



V ^(1) ^(1) J 



such that detF^^^ — detW^^^ , where the components of the matrix 



(1) 



(13) 
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are the Poisson brackets with respect to (g, tt) of the primary constraints . Hence, 
because / is regular, the regular or singular character of F^^^ is defined by the matrix 



(1) 



^ -R^,^^ 



(14) 



Thus, F^^j will be regular only if the Hessian R^^^ is regular, in which case the symplectic 
matrix can be constructed explicitly: 



1-1 



-f-'M,,,R-^ 



V 







R 



(1) 



(15) 



and we can define the generalized bracket 



A L-^(1)JaS 



1 ^^(Q^J 



^(1) 



(16) 



For dynamical functions that only depend on the (g, tt) variables it reduces to the Poisson 
bracket. On the other hand, if the Hessian R^_^^ is singular the matrix F^^^ is also singular 
and we cannot construct such a generalized bracket. To overcome this we must take into 
account other constraints, as we shall see in the following section. 

Here it is convenient to define the classification of the constraints we will use. We call 
primary constraints the ones discussed above, which are related to the transition to a first 
order theory. Secondary constraints directly arise from the zero modes of the Hessian of the 
original Lagrangian, and finally possible additional constraints are characterized as tertiary 
ones. 



III. THE SECONDARY LAGRANGIAN 

A non-regular Hessian R^^^ implies a singular F^^^, which is the signature of secondary 
constraints. These can be made explicit by the application of the F^^^ null eigenvectors on 
the equations of motion corresponding to the Lagrangian L^^^^{Q^^^). Thus, if F^^j has m <n 
zero modes i^^*^^, r = 1, ...,m, we have m secondary constraints: 

8 



%)(^,^>«>A,))-<,.^ = 0. (17) 



The zero modes are given by: 



V 



;., = (-<MJ-\u\-vl)^ (18) 



where the first component corresponds to the variables (g, tt), and and are the solutions 
of: 

<^a)=0 , = . (19) 

Introducing these vectors into Eq. (|1^ we get 



(2) o^ Qq. 0. Qq. Q^. ^g^.^ V ; 

and using the expressions (§) and (^ for H{q, vr, a) and fi*^^ respectively, we finally have: 

dL{q,a) d'^L{q,a 



(2) Oi 



dqt daidqj 



-a 



+ vJ'n' . (21) 
' i (1) \ > 



The last term may be disregarded because it is weakly null, and the remaining one is precisely 

an"- 

a Lagrangian constraint evaluated in g = a, which is independent of vTj, and thus -g;^ = 0. 
In the present schema the Lagrangian constraints of the original theory appear as secondary 
ones, a direct consequence of the singularity of the Hessian of this theory. 

Coming back to Eq. (^), and using the Poisson bracket defined by /"^ for the canonically 
conjugated variables (g,7r), the constraints n^"^^ may be written as: 

QJ' = v' {Q' , H} . (22) 

(2) oi L (1) ' J V y 

This shows that, although the primary constraints are considered in an analogous way by 
our procedure and the Dirac one, there is a relevant difference with respect to the secondary 
constraints. Our approach leads directly to a maximal set of linearly independent secondary 
constraints, in contrast with the Dirac method where a rather arbitrary linear combination 
of the constraints can be considered constraint. 
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Our following step is to incorporate these new constraints into the Lagrangian L^^^{Q^^^), 
using again a new set of velocity Lagrange multipliers, A* . This gives place to a new 

(2) 

Lagrangian 



=^(2)(Q(2)) =K)C) -^(Q(l)) ' 



(23) 



where Q = (g, vr, a, A , , A ). The symplectic matrix corresponding to this last Lagrangian 



is: 



(2) 



/ 





M, 



(2) 



R 



(2) 



T T 

, -M -R 

\ (2) (2) 



(24) 



with 



M. 



(2) 



f 9% \ 

dq dq 



R 



(2) 



da da 



(25) 



\ dn J 

Proceeding as in the previous step we construct the matrix W^^^, which has a structure 
analogous to W^^^ in Eq. (0)) but now with: 



(2) 



^ -^(2)^ 



(26) 



The components of uj^^^ are the Poisson brackets of the primary and secondary constraints, 



fijjj and with respect to {q,Tr). Taking into account Eqs. 

written: 



and (ESl) it can be 



(2) 



(2) 



(2) 



V 



^(1) 



do. 



(2) 



dq 



dq 





(27) 



In contrast with the primary Lagrangian stage, the matrix R^^^ is now rectangular. Our 
main problem is to study when F{2) is singular or regular, which is defined by ^E'(2)- This is 
considered in the two following sections. 
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IV. A REGULAR *(2) MATRIX: THE RESTRICTED PHASE SPACE AND THE 

GENERALIZED BRACKETS 



In this case F/o) is regular and hence F ^ exists. As the matrix ^ ^ exists, we can 

K^) O (2) (2) ' 



construct the matrix 



^\2) ^(2)^(2) ^(2) 



(2) 



(2) / 



(28) 



and easily obtain F^^^ by considering a block decomposition of J^- In fact, if we write 



(2) 



(^(2) )ll (^(2) )l2 ^ 



v 



(29) 



~^^(2)X2 (^(2))22 y 

where (\l/~^^)j^ and {^~^^)22 ^"^^ nxn and (n + m) x (n + m) antisymmetric matrices respec- 
tively, while (^~|^)i2 is a (n + m) X n rectangular matrix, the expression for F~^ becomes 



(2) 



(2) V (2) . 



(*r2)')-^(2)/"' 



(2)V (2)^12 



(2) V (2) . 
"(^(2) )l2 



V 



(30) 



1^(2) 22 -"'-'(2)^ V^(2)/12 1^(2) /'22 

Observe that in the block (1,1) of F'^ a structure quite similar to the Dirac bracket can be 

(2) 

recognized. 

This structure becomes more evident if u!(2) is regular, in which case one can perform an 
invertible transformation in order to write ^^^^ in a block diagonal form. This is supplied 
by the idempotent matrix 

\ 



S 



V 



/ 

-cu-^R" I , 

(2) (2) / 



(31) 



such that: 



(2) 



^ s s 

^ {2) — ^ ^(2)^ 



V 



uj 



(32) 



(2) / 
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where 7^^^ = R^^^u~^R^^y In this case 7^^) defines the regular or singular character of F(2). It 
is worth remarking that det 7^^) 7^ is equivalent to the condition of maximal rank for the 
matrix of the partial derivatives of the Lagrangian constraints with respect to the coordinates 
and the velocities, which warrants the soundness of the constraints 0. Consequently, in a 
well defined constrained Lagrangian theory 7^^) niust be regular. This allows us to obtain a 
more explicit expression for \I'^^^ and then construct out of the matrix R^^^ and 



( 



(2) 



where 



(2) (2) •' 



(2) 



-/-^M^^cj-^i?^ 7-1 

(2) (2) (2) '(2) 

7(-2) 
u-'R'' 7-1 

(2) (2) '(2) 



(2) 



(2) 



-/"'^,2)^(2) 
-7(2) ^(2) "(2) 



-'yjR,..ujJ 



C. 



(2) 



(33) 



(2) 



UJ 



(2) 



uj-^R"" -f-^R,,,u-^ , 

(2) (2) '(2) (2) (2) ' 



(34) 



T 

which satisfies the relation C,„,i? = 0, meaning that C,„, is constituted by the zero modes 

(2) (2) ' ° (2) -J 

T m 1 

of R^^^ . Therefore F~ has a structure analogous to the one that defines the symplectic form 
in the Dirac method, but extended to the whole configuration space, including the auxiliary 
variables. 

Thus, when uj(^2) is regular we can define a generalized bracket given by: 



dE{Q 



(2)' 



(2)^ 



dQ 



„i dG{Q 



(2)> 



(2) J AS 



(2) 



dQ 



(2) 



which in particular for the variables and Aj gives: 



(35) 



kl 



(36) 
(37) 



This case includes the systems with only second class constraints according to the Dirac 
method, and the phase space is a uniquely defined subspace of the primary one. 

It is also possible to have \l/(2) regular with u){2) singular. In this case we can redefine the 
constraints such that we can rewrite 
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, ,rxr 
^(2) 



( ~sxs n(' — s)xs 



Qsx{r—s) g(i — s)x{r—s) 



where r is the dimension and s the rank of a;(2) . Prom here on we can follow the preceeding 
procedure, but with U(2) instead of a;(2). This will lead to expressions similar to the ones 
already obtained in the case of a;(2) regular. 

V. A SINGULAR *(2) MATRIX 

When ^^2) is singular, F^^j is also singular, which makes it impossible to get a generalized 
bracket structure at this level. This may happen as a consequence of the existence of some 
still hidden constraints, which we call tertiary ones, or of a gauge symmetry of the secondary 
Lagrangian. In this case F^^) ^^^^ have the set of orthonormal zero modes {v^^^^ } 

v^^F^^^^O, s^l,...,m'<m, (38) 

where we use the symbol ~ to signify that this is a weak equation. The explicit form of v 



s 

(2) 



IS 



v;^ = {z^Ml^f-\f,z^), (39) 



where z'^ and y'^ satisfy 



z'Rl^O , z\^^+y^R^^^=0. (40) 



A. Tertiary constraints 

In an analogous way to secondary constraints, we can also have tertiary ones given by 

^:3)-<.)oag;;-4){^^)'^} + ^%)- (41) 

We have new constraints if H} are not linear combinations of the ^2^^^ and 0^2) ■ 

In this case we must iterate the procedure to include these constraints in a new Lagrangian 
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^(3)(Q(3))=A;<3,-W^(^a,)' (42) 
where Q^^^ = {q, q,a, \^^, \,^^, \.^^) . 

B. Gauge Symmetries 

When no new constraints arise from the above equation, the remaining zero modes give 
rise to the following symmetry transformations for the equations of motion (see Appendix 
C) 

^'Qi.,=K^^^ (43) 
which in terms of the coordinates of the extended configuration space is: 

6^Q = ez^Ml^r' = ezl^^{ni^,Q}, (44) 
5'ai = ey' , (45) 
6'X]=ez-^^. (46) 

This is a symmetry transformation of the secondary Lagrangian, but not necessarily a sym- 
metry transformation of the original one. One may wonder under which conditions the 
above coordinate transformations can be promoted to symmetries of the original action. To 
answer this question we can consider the variation of the action induced by transformation 
(1^ ), taking into account that 

zi^{n^,),H} = ui,^n^,) . (47) 

Thus we finally get 

rt2 r /an \i 

(48) 



r dt€ 



The surface term in the case of internal symmetries will cancel only if the constraints are 
first order homogeneous functions of the canonical momenta tt, as happens in the Dirac 
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method. It is worth remarking that it is possible to add a term to the time derivative of 
in order to cancel the residual terms, which are all of them proportional to the constraints 

d,e - Dte = d,e - z^^ (f/^ fi^^, + y^fi^,,) e . (49) 

Here we have assumed that the zero modes of are orthonormalized and there is no sum 

(2) 

on s. In this way, by introducing this covariant derivative weakly equivalent to the ordinary 
one, we obtain a symmetry holding on the whole of the configuration space, whenever Eq.(p8D 
is strongly satisfied. 

The relations (^) may be fulfilled in the following two situations. 

a) 2;"* = and y'^ is a zero mode of R^^y Given that R^^^ is singular, this implies that 
the matrix {—^) has a zero mode. In this case v^^^^ = {0,y^,0), and does not lead to new 
constraints but generates the transformation symmetry 

6Q = (50) 
Sa = esV' (51) 
^K, = (52) 

under which the Lagrangian variation is: 

OVl OH 

= ^"^9^^ + '^"^a^ = + ^(^'^^ = • ^^^^ 

h) and so there are non-trivial zero modes of i?^^. If tu^^) is regular and 7^^) is 

singular, the relations (^ ) become: 

= y'^i^^^;.' ' = ' (54) 

so that the zero modes can be written 



V 

(2) 



(^^^(2)^(■;;<)/"^^^-^^^(2)^(■;;)- (55) 

and the associated constraints are 

%,=^t%; + v'R,.,^;,'{^,H}- (56) 
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The first term is a linear combination of tlie primary constraints, and tlie second one contains 
a matrix whose components are the Poisson brackets {Q^^^H}, a hnear combination of 
the secondary constraints, and {Q^^^H}. The primary and secondary constraints have 
already been considered, and thus the tertiary constraints will only correspond to a given 
combination of the Poisson brackets {Q^^^H}. In general, in this case we must carry on 
and include these new constraints in a tertiary Lagrangian, by iterating the preceeding 
construction. Tertiary constraints may arise only when the set of zero modes of 7^^, is 
bigger than the set of zero modes of R^^^ , and whenever the time evolution of the secondary 
constraints, {Q^^^,H}, does not reduce to a linear combination of fi^^^ and Q^^y If no 
new constraint arises, the transformation generated by (|^) leaves the equation of motion 
invariant on the constrained subspace. 

In the simplest situation where both sets of zero modes are equal, we have v^R^^^ = 
and recover the situation described in the preceeding subsection. 

VI. SOME EXAMPLES 

To clarify the use and possibilities of our method, we consider here two simple examples 
which have been widely discussed in the literature. Both present pathological features from 
the point of view of the Dirac approach. The first one is a regular Lagrangian, but with 
terms of higher degree in the derivatives. This causes a bifurcation in the canonical structure 
according to the usual approach, which gives place to a problem with multivaluated boundary 
conditions. Our approach allows us to construct a symplectic formalism with physical sense, 
which also gives a consistent quantum theory. The second example contains a first class 
constraint, according to the Dirac classification, which has a weak null bracket with every 
dynamical function, and thus does not admit a gauge fixing. In this case our approach 
permits a deeper understanding of the origin of the pathology. 
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A. 1.- Lagrangian of higher order time derivatives 



Let us now consider an example with higher order time derivatives, which has been 



studied in Ref. fTO], defined by the Lagrangian 



m = \{qY~\k{q? (57) 
The equation of motion in the configuration space is: 

(2>q^ -k)q = 

which imphes q = cte., i.e. it describes a one-dimensional free-particle. 

When a Legendre transformation is used to construct a canonical formalism, this La- 
grangian leads to velocities which are multivaluated functions of the canonical momenta, 
making the classical motion unpredictable since at any time one can jump from one branch 
of the Hamiltonian to another. We now analyze this system in the light of our procedure. 

Proceeding as proposed in the first section, we introduce the auxiliary variables a and p 
to obtain a first order time derivative Lagrangian i^(o), 



L(o)(g,p,a) =p{q-a) + ^a^-^ka^ . (58) 



Its equations of motion reveal the constraint 

n = ^ = a'- ka-p. (59) 
da 

By incorporating it to the Lagrangian we get 

{q, p, a, X) = pq + {a^ - k a - p) X -H{p, a) , (60) 
where the Hamiltonian is 

H{p, a) = — ^ + ^ka'^ + pa . (61) 

Arranging the coordinates in a four-dimensional vector = {r = q — X,p,a, X), one easily 
obtains the symplectic matrix F^^^ 
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(1) 



0-10 
10 
Sa^-A; 
k-3a^ 



Its inverse defines a generalized Poisson bracket 



^1 ^^^(Q,.) 



so that the equations of motion can be written as 



Q,, = {Qt,Hr. 



(62) 



(63) 



(64) 



They lead to the classical solution q{t) = cte, consistent with the original Lagrangian dy- 
namics, but in fact the configuration space is now two-dimensional, without constraints, 
and therefore we have two degrees of freedom and we need four boundary conditions. This 
establishes a difference with the usual description for a free particle, given by a Lagrangian 
quadratic in the velocity. The equations of motion (B3) can be displayed as: 



(65) 



r = a 

\ = — ka — p 
a = 
p = 

which implies that the four phase-space coordinates are constants of motion, and hence Q 
= f + A. is also a constant of motion. But now it is not sufficient to give, for example, the 
position and the velocity at a given time: 



qo = r + X , 

qo = + {1 - k) a - p , 



(66) 



because there is not enough information to determine a point in the phase-space. 



There is another remark that should be made. The phase space coordinates Qj^^ are 
not canonical, in the sense that the brackets of the fundamental variables are not equal to 
one or to zero. If we try to obtain a canonical structure through a change of coordinates 
it is necessary to consider a non-linear transformation, such as ={r = q — X,p,a = 
a^—ka,X). These new coordinates satisfy canonical brackets, but this type of transformation 
also introduces bifurcations in the phase space with a correlated multivaluated Hamiltonian. 

Although the Qj^^ coordinates are not canonical, the dynamics described in their terms 
has a symplectic structure. This is given by Eqs. (|62| , |63|) , which define the evolution of 
the system through the generalized Hamiltonian (^) and the equations of motion (|64|). 
This symplectic structure can be used to construct a quantum theory by using the canon- 
ical procedure given by {A^B}* —i[A,B], where A and B are the quantum operators 
corresponding to the classical functions A and B. This procedure leads to: 

[f,p] = z, (67) 
[a, A] = i(k-?,a^y^ , (68) 



with all the other fundamental commutators null. A representation for the fundamental 
operators that satisfies these commutation relations is: 

f = r , d = a , (69) 

p = -i-^ , X = -i(k- 3a^) ^ ^ . 
or ^ ' oa 

The time derivatives of the fundamental phase space operators are given by: 

r = —i[f,H] = a , 

d = -z[a,^]=0, (70) 
P = -i[P,H] = , 

A = -i[X, H] = - (k- 3a^) - + ka) = -a + ^ 



Sa^ — k 



which implies that p, a, ^ and ^ are constants of motion. For the coordinates of the config- 
uration space, r and A, and their velocities we have the following commutation relationships: 
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^ dr 

T 

d\ 



0, 



- 



-1 



d_ 

da 



—a + 



Sa^ — k 



k — 



6a 



p — a (a^ — fc) 
{3a^-kf 



(71) 



(72) 



In particular the velocity of the particle 



dq df ^ dX p + 1o? 
dt dt dt 3q;^ — k 



(73) 



is a constant of motion that can take any real value, as is characteristic for a free particle. 
But the commutator of the position and the velocity is here: 



^ dq 









r + A, 



p + 2a^ 
3a^ — k 



3a^ — k 



1 + 



Qa 



+ a {ha^ 



k 



+ ^ ■ (74) 



(3a2 - ky 

This implies that the uncertainty relationship between the position and the velocity is dif- 
ferent from the one resulting when the Lagrangian is quadratic in the velocity. This result 
clearly shows that we have a new quantum theory, but with the same classical limit than 
the quadratic one. 



B. 2.- Lagrangian with a pathological Dirac constraint 

The following Lagrangian was proposed by Cawley: 

^ = q\q2 + gag? • (75) 

The construction of a canonical formalism for this system is cumbersome. This and other 
related examples have motivated several discussions and ad hoc modifications of the Dirac 
method. The Lagrangian equations of motion give: 

= arbitrary constant , 
52 = , (76) 
§3 = arbitrary function of t . 

If we try to apply the Dirac method we find a primary constraint, = 0, and a secondary 
one, qf = 0. This latter has been called ineffective [0 or a fourth class |§] constraint, 
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for which it is impossible to construct a gauge fixing constraint. This prevents us from 
determining a physical subspace of the phase space with a symplectic structure. To overcome 
this difficulty we could try to use a linearized constraint, ^2 = 0, instead of the secondary 
one. Prom an algebraic point of view both constraints are equivalent, but this last one leads 
to an additional constraint, pi — 0. Thus we would now have three first class constraints, 
7^3 = 0, ^2 = and pi = 0, for which it would be necessary to introduce a gauge fixing. Each 
gauge fixing would define a configuration characterized by: 



and so it would not reproduce the whole solution space of the Lagrangian equations. For 
this reason this Lagrangian is considered a pathological one from the point of view of the 
Dirac method. 

Now we will apply our method to this Lagrangian. To simplify the discussion we use the 
shortcut provided by the algebraic elimination of cyclic variables. The first order Lagrangian 
is: 



qi — a given constant , 



q2 = 



(77) 



qz — a given function of t , 



3 



(78) 



i=l 



with 



3 



T^iUi - aia2 - qsql ■ 



(79) 



1=1 



The equations of motion of the a variables give the constraints: 



ni 





(1) 



^-^(0) 



= -TTs . 



= «! - 7r2 , 



= q;2 - TTl , 



(80) 
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We introduce these constraints into the Lagrangian by means of a set of velocity Lagrange 
multiphers A, and thus we obtain the primary Lagrangian L(i) 

3 3 

^(1) = E 9^ - E ^(1) . (81) 

i=l i=l 

The Hessian has a zero mode leading to the secondary constraint: 

n,:, ^ ^ ^ . (82) 
Again, we incorporate this constraint to the Lagrangian as before: 

3 3 

= E ^» ~ E ^(1) '^i-QiV-H , (83) 

1=1 i=l 

and now the rectangular matrix i?(2) is 

R{2) = R(i) I Osxi . (84) 

It has two zero modes satisfying the Eqs. (|^). One of them is basically the same as in the 
previous step, = (0, 0, 1, 0) , and gives rise again to the constraint Q(^2) 

Q},) = /f} = ^ = -ql = ^ . (85) 

The other zero mode, z"^ = (0, 0, 0, 1), yields the tertiary constraint 



which is weakly zero by virtue of fi(2)- Hence no new constraints actually arise because 
of these last two zero modes and we may conclude that these zero modes give rise to a 
symmetry of the equations of motion. 

In fact, the symmetries just found correspond to the last set of equations of motion, 
which include all the constraints, and do not necessarily correspond to symmetries of the 
original Lagrangian. The question is now whether it is possible to extend them to the 
original Lagrangian or not. At this point the symmetry generated by the second zero mode 
poses a problem. Observe that if we try to carry the expression of fi^g-, into the form of Eq. 
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(PI), the coefficients U^,^-^, given by Eq. p7|), become ill defined on the constrained subspace, 
thus preventing the extension of this symmetry to the whole configuration space, as one can 
see through the definition of the covariant derivative given by Eq. (^Q]). 

In summary, the symmetry generated by the ffist zero mode can be extended to the 
whole phase space by constructing an adequate covariant derivative, which enables us to 
introduce a gauge fixing without altering the original Lagrangian equations of motion. But 
this is not the case for the symmetry generated by the second zero mode, which is strictly 
valid only in the constrained subspace. Any additional condition introduced to make the 
presimplectic matrix regular only corresponds in fact to a gauge fixing in this subspace. 
Therefore, this additional condition will maintain the equations of motion in the subspace, 
but will alter them in the whole space, and leads to a set of equations of motion different to 
the one given by the original Lagrangian. 



VII. FINAL REMARKS 

In this paper we develop a general theory of singular Lagrangians and their constraints. 
It leads to a systematic and purely Lagrangian approach for constructing a symplectic struc- 
ture, thus avoiding the use of a not always well defined Legendre transformation and the rise 
of many associated pathologies. Our approach is based on an enlargement of the original 
configuration space (the q space), which makes use of two sets of auxiliary variables. One 
set allows us to rewrite the original theory as a ffist order one (the {a, n) variables) and 
the other contains the Legendre multipliers (the A variables), which we use to incorporate 
all the constraints at the level of Lagrangian equations of motion . These last variables are 
introduced as velocities in the enlarged Lagrangian, such that they acquire the role of con- 
jugate coordinates of the constraints. In this way we have a well defined first order regular 
dynamics in an extended configuration space, which warrants that the q coordinates satisfy 
exactly the same equations of motion than in the original Lagrangian. This procedure al- 
lows us to disentangle all the mechanisms hidden in the Dirac-Bergmann approach, and to 
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identify the source of the pathologies that appear there. 

Our construction gives place to three kinds of constraints. In all the cases we have a 
set of primary constraints, which resemble the definition of the canonical momenta. They 
characterize the enlargement of the configuration space necessary to pass from the original 
theory to a first order one. For a regular Lagrangian they are the only constraints that 
appear and our construction stops here. When the original Lagrangian is singular a new 
set of constraints appear, which we call secondary constraints. They are associated to 
the zero modes of the Hessian of the original theory. Finally, there may be additional 
hidden constraints, not directly related to the zero modes of the Hessian, which we have 
characterized as tertiary ones. 

These constraints can hold in the whole configuration space or in a constrained subspace. 
According to this they can be classified in two categories. One of them corresponds to what 
can be called strong constraints, i.e. the set of independent constraints that act on the whole 
configuration space and give place to a reduction of the extended space. They contain the 
second class constraints of Dirac. 

The other category is constituted by the weak constraints, i.e. constraints that hold 
only on a constrained subspace. Strictly speaking, the latter are the generators of symmetry 
transformations in such a constrained subspace. Here there are two possibilities, according 
to whether they close an algebra or not. In the first case the symmetry can be promoted 
to be valid on the whole configuration space by defining a covariant time derivative for the 
parameter of the transformation, weakly equivalent to the ordinary one. Hence we can give 
a gauge fixing that does not alter the equations of motion of the original coordinates g, and 
we can construct a symplectic form. In the second case, when the constraints do not close 
an algebra, the symmetry holds only on the constrained subspace and cannot be extended to 
the whole space. In this case we can give a gauge fixing in the constrained subspace, but it 
will alter the equations of motion of the original coordinates and it thus becomes impossible 
to construct a symplectic form without altering the original equations of motion. 

In summary, when there are strong constraints or weak ones that close an algebra our 
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approach allows us to construct a consistent symplectic formalism. Our brackets are defined 
on the whole extended space, and correspond to the Dirac brackets in the {q, tt) subspace 
when they exist. Furthermore, we can perform this construction even when there is no 
well defined Hamiltonian, as in the example discussed in Section IV-1, which renders our 
approach more general than the Dirac method. 

When there are weak constraints that do not close an algebra, we cannot construct 
a meaningful symplectic formalism consistent with the original equations of motion. The 
Lagrangians that lead to this case are pathological within the Dirac method, and also in 
other approaches. 

In contrast with the Dirac method, in our approach the whole set of constraints is 
unambiguously defined and there is a clear connection between the constraints and the 
symmetries of the Lagrangian. Some of the main contributions of our work are that it gives 
a simple and unified framework for understanding the origin of the pathologies that appear 
in the Dirac-Bergmann formalism, which are explained by the impossibility to promote 
a given symmetry from the constrained subspace to the whole configuration space, and 
that it brings a more general approach for constructing a symplectic formalism, even when 
there is no Hamiltonian at all in a canonical sense. This last point opens the possibility of 
overpassing the limitations of the canonical quantization, and performing an algebraically 
consistent quantization on the basis of this non-canonical symplectic algebra. 
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APPENDIX A: THE CANONICAL HAMILTONIAN 

In the very particular case in which all the cyclic variables a, can be algebraically elim- 
inated, we are finally faced with a Lagrangian of the form: 
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(0) 

with the equation of motion given by 

jabQb ^ 

where = {qi,'^i) and 



(Al) 



(A2) 



(A3) 



a symplcctic matrix in the (q, n) space. This situation corresponds to a regular Lagrangian, 
and the equations of motion can be written 



A{Q) = {A{Q),H{Q)} 



where H{Q) is the usual canonical Hamiltonian, in terms of the Poisson brackets 

iF(n) rron-^^(^) f-^ ^^^^^ 



(A4) 



(A5) 



APPENDIX B: THE DIRAC HAMILTONIAN 



Let us consider a Lagrangian bilinear in the velocities: 



L = -mijqiqj - V{q) 



i — 1, n 



(Bl) 



where the matrix niij is singular. The procedure presented in Section II produces the first 
order Lagrangian: 



-^(0) = ^iQi — T^iC^i + H , 



(B2) 



with 



which leads to the following primary constraints: 
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(B3) 



^(1) = rUijCij - TTj . (B4) 

We assume that the range of rriij is r < n, and introduce the operator U that diagonahzes 
the mass matrix 

iU^mU)j = XAj , (B5) 
and transforms the coordinates according to 

a a^U'^aU, (B6) 

TT ^ TT = U'^TtU . 

After this transformation the constraints become 

TTs = Xs&s s = 1, ...,r (B7) 

TTp = p^r + l,...,n 

We can use these constraints to ehminate the first r cychc variables ag, but the last n — r 
ones ap will remain arbitrary. In such a way the Hamiltonian becomes: 

H = ^Xj^TTsTTs + V{q) + apfcp , (B8) 

which in terms of the original variables (q's, tt^) and the arbitrary functions ap can be written 

H = ^ ''' TTsTTt + V{q) + apTTp , 

which is the Dirac Hamiltonian. 

APPENDIX C: CONSTRAINTS AND SYMMETRIES 

Given a first order Lagrangian: 

L{Q,Q)^f{Q)Q-H{Q), (CI) 

the form of the corresponding equations of motion is: 
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FiQ)Q-{V^H){Q) = 0. (C2) 

If the matrix F{Q) has zero modes v^, v^F{Q) = 0, from the equations of motion we obtain 
an associated set of constraints 

n^{Q)=v^{V^H){Q). (C3) 

The transformation of coordinates generated by these zero modes 

6Q = v'er , (C4) 

where are time-dependent parameters, produces the Lagrangian variation 

6L = ^6Q = erv'{F{Q)Q - {V^H){Q)) = erOT . (C5) 
oQ 

The relation (|C5|) shows that the Lagrangian is weakly invariant under displacements 
in directions orthogonal to the gradient of the potential. In a Lagrangian formalism the 
constraints are part of the equations of motion. In this sense this sjmametry is valid on a 
subspace determined by the solutions of the equations of motion, and can be called a "weak 
symmetry", in contrast with the symmetries independent of the solutions of the equations 
of motion which can be called " strong symmetries " . When the constraints close an algebra, 
this symmetry is related to a strong one, obtained by a simple redefinition of the derivatives 
of the parameter e. 
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